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The influence of a constant uniform magnetic field on the thermodynamic properties of a partially 
ionized hydrogen plasma is studied. Using the method of Green's function various interaction 
contributions to the thermodynamic functions are calculated. The equation of state of a quantum 
magnetized plasma is presented within the framework of a low density expansion up to the order 
e 4 n 2 and, additionally, including ladder type contributions via the bound states in the case of strong 
magnetic fields (2.35 x 10 5 T « B < 2.35 x 10 9 T). We show that for high densities (n w 10 27_30 m~ 3 ) 
and temperatures T « 10 5 — lCf'K typical for the surface of neutron stars nonideality effects as ,e.g., 
Debye screening must be taken into account. 
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O I. INTRODUCTION 

The calculation of the equation of state (EOS) of a multi-component quantum plasma consisting of charged particles 
interacting via the Coloumb potential is of theoretical interest as well as of practical relevance, e.g., for astrophysical 
systems such as stars. The aim of this paper is to derive a low-density expansion for the equation of state (EOS) 
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of a two-component plasma embedded in an external constant magnetic field. This problem was recently tackled by 



Cornu [jj and Boose&Perez Q who derived a formally exact virial expansion of the EOS by using a formalism which 



, ^ \ is based on the Feynman-Kac path-integral representation of the grand-canonical potential. 

In this paper we will employ the method of Green's function. As the calculations are carried out for a nonrelativistic 
quantum system, we restrict ourselves to magnetic field strengths B < B re i which is given by B re i = m\c 2 / (eh) rs 
4.4 x 10 9 T. Further we will use an expansion of the magnetized plasma pressure in terms of the fugacity z — e^** 
to obtain the EOS of a weakly coupled magnetized plasma. Thus we can derive explicit expressions for various 
contributions to the quantum second virial coefficient. Though the formalism is formally valid only for low densities 
the obtained explicit expressions are appropriate even at sufficient high densities as the magnetic field increases the 
domain of classical behavior towards higher densities. The second virial coefficient contains both scattering and bound 
state contributions of two-particle states. Being interested in the thermodynamic properties of quantum magnetized 
plasmas the influence of the magnetic field on the energy eigen states of a two-particle states has to be taken into 
t*"*- ■ account. 

Usually the magnetic field is measured by the dimensionless parameter 7 — Tluj c / (2Ry) = B/Bo where hu> c is the 
electron cyclotron energy, Bq w 2.35 x 10 5 T, and Ry — e 2 /(87reoas) w 13.605eU is the ionization energy of the 
field- free hydrogen atom. Whenever 7 > 1, i.e. the cyclotron energy is larger than the typical Coulomb energy, the 



, structure of the hydrogen atom is dramatically changed. This problem has been approached by several authors 



Using the results of these authors we study the influence of bound and scattering states on thermodynamic properties 



of magnetized plasmas. 



Recently the problem of ionization equilibrium of hydrogen atoms in superstrong magnetic fields (7 » 1) was 
considered by Lai&Salpeter ||. They proposed an ideal Saha equation of a hydrogen gas including bound states but 
neglecting screening effects and scattering contributions to the second virial coefficient. Using the EOS obtained in 
our derivation we construct a modified Saha equation which takes into account nonideality effects as well. 

The paper is organized as follows. In section II, we discuss the method which is used to calculate thermodynamic 
functions and derive analytical results for the scattering contribution in section III. An approximate result for the 
C$ I bound state contributions is given in section IV and the equation of state is presented in section V. Finally, we use 
our results to derive a generalized Saha equation and compare the degree of ionization with the results of the ideal 
Saha equation in section VI. 



II. FUGACITY EXPANSIONS OF THE THERMODYNAMIC FUNCTIONS 

We consider a two-component charge-symmetrical system of N spin half particles of charge (-e) and mass m e and 
N spin half particles of charge e and mass m,. In general, the total pressure can be split into ideal contributions and 
interaction contributions 
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P = Pid + Pint ■ 



(1) 



The pressure and the particle density of an ideal plasma in a homogeneous magnetic field B = (0, 0, Bq) are given by 
a sum of Fermi integrals over all Landau levels n 



^ = fcT E |f E'/fM > ^ = £lf E'-Mino ( 2 ) 

=0 a a n=0 



(x a — hujc I (2kT) with = \e a \Bo/m a , A a = h/ \j2-Km a kT, and z° = exp — n%U)^)}). The prime indicates the 
double summation due to the spin degeneracy except for the n — level. 

The interaction part of the pressure for sufficiently strong decaying potentials may be written in terms of a fugacity 
expansion 

Pijp- Pid) = ^ZaZbBab + ^ Z a Z b Z c B abc + ... , (3) 
ab abc 

where we have introduced the modified fugacities 

~ 2 x a 



Af L tanh(xa) 



(4) 



In the limit of small densities we have z a — > n a . We focus on the calculation of the second virial coefficient B ab which 
is defined by 

Bab = J_ ( A a tanh(a; a ) \ / tanh(a: fe ) \ ^ ^-0^ _ , (5) 



2n V 2 x a J V 2 Xfo 
H^ b is the Hamiltonoperator of the two particle system with the interaction potential V^&(r) 

H X ab - ( iPa ~ €aAa)2 + li a B B Q a z \ + ( (Pb ~ ebA " )2 + n b B B a z ) + XV ab (r) , a z = -1, +1 (6) 



2m a J \ 2m b 

■\ - of the noninteracting system. The additive term fi B B a z describes the coupling between the intrinsic 
magnetic moment (pL a B = e a h/ (2m a )) of the charged particles and the magnetic field. However, in the case of particles 
interacting via the Coulomb potential V a b(r) = e a ef,/(47re | r ra — r& |) the second virial coefficient defined by Eqs. ([]]](]) 
is divergent. In order to obtain a convergent expression one has to perform a screening procedure. Such a technique 
is well established in the zero magnetic field case and can be easily extended to the nonzero magnetic field case. 
This program was also carried out by Cornu Q and Boose&Perez || who used the Feynman-Kac formalism to derive 
a virial expansion for a magnetized multi-component system. Using the methods as described in the convergent 
second virial coefficient of a plasma may be split into a scattering and bound state contribution. In contrast to the 
zero magnetic field case an exact calculation of the convergent second virial coefficient in terms of scattering phase 
shifts is very complicated. Therefore we will give a perturbation expansion of the scattering part in terms of the 
interaction parameter e 2 up to the order e 4 and use an approximate expression for the bound state part which is valid 
in the case of strong magnetic fields (7 > 100). We may employ the method of Green's function. The starting point is 
the observation that the equation of state is connected to the average interaction energy (XV a b) by a charging process 

1 f 1 dX nr x 

P-Pid = -Q J —{Vab)\, (7) 

fi is the volume of the system. Taking into account many body effects thermodynamic functions may be expressed by 
a screened potential V^ b . By this method the divergencies due to the long range Coulomb force are removed. Then 
the pressure is given by the equation 

1 p ^ a \ r 

HP- Pid) = go £ / — / dl d2 (V ab (12X) G a (ll) G b (22) + V° b (12\) n ab (121++2+A)) . (8) 

ab 

Here the first term is the Hartree approximation given in terms of the free particle Green's function G a (ll) and H a b 
denotes the polarization function. For low density systems it is necessary to calculate bound state contributions to 
the thermodynamic functions. Therefore we apply the ladder approximation for H ab 
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(9) 

To avoid double counting we have introduced the operator P3 which subtracts contributions of the order V^ b and 
(VA) 2 . We may divide p^t into a bound state contribution pf°" and a scattering state contribution pf%% u 

n . — bound , scatt /-,r,\ 
Ant — Pint 'Pint • \ lv I 

In the case of a Coloumb potential this division is not trivial as the atomic partition function is divergent due to the 
infinite number of bound states at the continuum boundary. This problem has been extensively discussed in the zero 
magnetic field case Q . One can solve this problem in a natural way by introducing a renormalized sum of bound 
states 

bound _ ~ ~ p rtbound fill 
Pint — z eZl*3 D ab , [LL) 

where at zero magnetic field B h ° b und is given by the Planck-Larkin partition function . This division is somewhat 
arbitrary but guarantees the convergence of the bound state partition function even at vanishing magnetic field. We 
mention that this division does not affect the results of the thermodynamic potentials. 



III. SCATTERING STATE CONTRIBUTION 



We consider all diagramms up to the order e in the interaction parameter. A diagrammatic representation of the 
perturbation expansion takes the form 




These diagrams arc the Hartrcc term, the Montroll-Ward term, the Hartree-Fock term, and the exchange e 4 term, 
resp ectively. The solid lines represent the uncorrelated Green's function for a charged particle in a magnetic field 
|pl| . Hence our calculations are valid at arbitrary magnetic field strength. The divergence of the Montroll-Ward 
graph is avoided by introducing a screened potential line. The screened interaction potential V s is evaluated in the 
random phase approximation V s (q, oj) — V(q)/(l — V(q) H RPA (q, lu)). At low densities V s can be approximated by a 
statically screened potential V s = e 2 /(e [q 2 + n ,}) with n 2 = (e 2 /e ) U RPA (0, 0) = (3 (e 2 /e ) (z e + h)- In the following 
calculations all results are obtained by setting the distribution function fo(u>) = e' 3M e _/3aj , i.e., in the nondegenerate 
limit n\ 3 tanh(x)/a; -C 1. The Hartree term vanishes due to the electroneutrality. 



A. GREEN'S FUNCTION FOR THE MAGNETIC FIELD PROBLEM 



In this section we represent the uncorrelated Green's function for a charged particle moving in a constant magnetic 
field in a closed form. The Green's function is the solution of the equation of motion (using symmetric gauge and 
setting h = 1): 

" lr {x2 + y2) + T Lz ~ > 1bB(Jz + *w) G ' (R ' T) = m 5{T) ■ (13) 

G'(r, r',T) can be expressed in terms of the correlation functions by G'(r, r',T) = 0{T) G>(r, r', T) + 
0(— T) G<(r, r',T). The prime denotes the particular choice of the gauge. Both G> and G' K satisfy the homoge- 
neous counterpart of Eq. (|l3|) . According to Horing Jll| , for arbitrarily chosen gauge they can be written as, 

G { > } (r,r',T) = J g | ^jj"^ | exp(-zcT) J°° dT> exp(iuT') A(r, r', T') (14) 

with 
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A(r,r',r') = C(r,r') ^ ^ exp( ?P R) exp 



—i hbB(7 z 



2m 



r 



r ex P 



.pI+pI 



tan —T 



(15) 



The gauge dependence of the Green's function is explicitly given in the factor C(r, r'). Noting that C(r,r') is only a 
function of R = r r' and that it obeys the relation C(r, r')C(r', r) = 1, this factor can be left aside in the following 
calculations. 



B. Hartree-Fock term 

First we calculate the Hartree-Fock term, which can be written in space time representation as 
/W = -^£ J -Tr (a) J^dld2V{l2)Gl{12)Gl{21+)5 ab . 



(16) 



a b 



The free particle Green's function G£(12) must now be replaced by Eq. (|l4|) . In the resulting expression all integrals 
can be computed exactly. The detailed calculation is given in Appendix A. Defining t; ab = e a e b / (4TreokT\ ab ) and 
Kit = h/ \/2m ab kT, m ab being the effectiv mass, we obtain the result 

0PHF ^^2jZa 2 X 3 a a S,aafl(Xa), (17) 



where we have introduced 



tanh(a; a ) cosh(2a; Q ) arctanhy'l 



tanh(x a ) 



X a COSh 2 (x a ) tanh(x Q ) 



(18) 



C. Montroll-Ward term 



Next we investigate the direct term of order e 4 given by the following expression 

ab 



PPmw = ^J2 fy^^) / dld2d3d4V° b (12)V ab (M)G:(23)G:(32)G° b '(U)G° b '(41). 



(19) 



Again, a detailed calculation may be found in Appendix B. Retaining only contributions of order z 2 we obtain the 
result 



3 — 

Ppmw = ~rz y~] —r l aZ b \ 3 ab £ b f 2 (xa, x b ) 



7T2 

~4 



(20) 



where f%{x a , x b ) may be written as, 



1 4 



2 7T 



h{x a ,X b ) = - + - / dty/t(l-t) (Va + Vb) 



arctanh^/1 - (y a + y b ) 
V 1 - {Vo. +Vb) 



(21) 



with y a ^ b = \ 2 aabb sinh(x a ,fci) sinh(a; a:6 (l - t))/(X 2 ab t(l - t)2x a ^ b sinh(a; a:6 )). The first term in Eq.@ is the Debye 
limiting law, while the second term gives a quantum correction. According to the Bohr-van-Leeuwen theorem the 
classical Debye law is not influenced by a magnetic field. 
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D. second order exchange term 



The exchange term of order e 4 is given by 

" x dX 



1 r r1\ (' 

2n^2j T Tr ( ff ) J dld2dZd±V ab (lZ)V ab {2±)G a (\2)G°(2Z)G a {M)G a {±l)5 ab . (22) 



0P, 

The result can be written in the form (Appendix C) 



0P+ = " E Z " A - & , (23) 



where fs{x a ) is given by an integral representation (|C4j) and can only be evaluated numerically. Therefore we propose 
the following fit expression for fs(x a ) 



_ cosh(2x a ) / tanh(cx a ) \ d arctanh-y/l - ^^fy^ 

COSh 2 (x a ) \ ( CX a) J !\ tanh(c"^~y 

V ( ra ») 

with the fitting parameters c= 0.8349 and d=0.9169. 

Finally, we may sum up all contributions up to the order z 2 e 4 . Collecting the obtained results (17 2C| , |23] ), the 
scattering states contribution to the pressure in this approximation may be written as, 



12tt 

ab 



where we have defined B^l att by 



BT tt = \^a b lxlU ab f 1 (x a )~^Xl b e ab f2(x a ,x b )~6 ab ^ln(2)\l b e ab h(x a )J . (26) 

The influence of these states on the thermodynamics will be studied in section V and VI. Finally, we note that this 
equation gives in the limit x a — > the exact zero magnetic field results (see ||). 



IV. BOUND STATE CONTRIBUTION 

According to Eq.(]|) we have for the bound state contribution 

Ppt° n r d = z eZl P z ^' PEm - (27) 

m 

where E m are the eigenvalues of H^ 1 . In Eq.(^) all terms up to the order e 4 with respect to the interaction 
parameter must be omitted. In order to calculate p h °^ nd the precise knowledge of the binding energies is essential. 
Therefore we briefly review the energy spectrum of the bound states and specify the approximations used in this 
paper. In contrast to the field-free hydrogen atom there is no exact solution for the nonrelativistic hydrogen atom 
at abritrary magnetic field strength. We focus on the astrophysical interesting strong field regime 7 > 1. Here we 
essentially follow the work of Lai&Salpeter H . 

The two-body problem has been investigated in the pseudomomentum approach |^|^|. The pseudomomentum 
K = X) Q (Pa ~ e a A a + e a B x r a ) is a constant of motion. Therefore one can construct a wave function with a 
well-defined value of K by 

tp(R, r) = exp[i(K + (1/2)B x r)R]<?!>(r) , (28) 

with the centre of mass-coordinates R = (miri + m2£-i)l(m\ + 1112) and the relative coordinates r = r-2 — ri. Then 
the Hamiltonian of the Schrodinger equation H<f){r) — (Hi + H2)<f>(r) = E nml/ KzK ± 4>(r) can be written in the form 
(setting A = 1/2 (B x r) and M = m e + rrii) 
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Hi = ^ + ^— (Bxrf+(- -) |B(r x p) - , (29) 

2m e i 8m e i \m e rrii J 2 iireor 

( m e \ hiu e r K 2 K 2 e , 

In this approach the spectrum is characterized by the Landau quantum number n of the electron, the magnetic 
quantum number m, the number of nodes v of the z wave function, and the pseudomomcntum K. In case 7 > 1 we 
can restrict ourselves to n=0. The energy eigenvalues read as [|| 

m K 2 K 2 

E 0mvKx K x = E mv + mHul— + 77V7 + (3!) 

rrii 2M 2M± 

E mv is the energy of a bound electron moving in a fixed Coulomb potential. For v = Q the states are tightly bound 
with binding energies 



E m0 = -0.32 ^ In 2 [Tr^—^t ) RV > (32) 
m e \ 2m + 1 m~ 



while for v > 1 the states are hydrogen-like and the eigenvalues are well approximated by 

E mv =-4 — Ry ^ = 1,2,3,4... (33) 
v{ m e 

for the odd states (i.e. v — 2v\ — 1) and for the even states (i.e. v — 2v\). The second term in Eq.(|3l|) describes a 
Landau excitation of the proton which is coupled to the electron quantum number m due to the conservation of total 
pseudomomentum. The atom can freely move along the magnetic field direction contributing the term K 2 /2M to 
the energy. Contrary to that the transverse motion is coupled to the internal motion by the term (e/M)(K x B) r. 
For magnetic field strengths considered here energy corrections due to this term can be computed by pertubation 
expansion with respect to the eigenstates of H\. Lai&Salpeter proposed an effective mass M± approximation of the 
transverse moving atom with 

Mi = M 1 1 + t J — — I , t » 2.8 , (34) 

^ 0.32^1n( 7 )y 

which we will use for simplification for all m-states. This energy correction is only valid for small pseudomomentum 
K± < K ±c where K ±c is defined by h 2 K\J(2M) « (0.32(M/m e ) ln(7)/(*7)) Ry but serves as a fair approximation 
for magnetic fields strengths B < 2.35 x 10 9 T. We note that due to the coupling of the intrinsic magnetic moment 
of the proton with the magnetic field an additional factor of (l + e~ 2:Ei ) arises in the bound state partition function. 
On the other hand, at magnetic fields 7> 1 and temperatures T sw 10 5-6 K spin excitations of the electrons can be 
neglected. 

Given the energy eigenvalues we can define a convergent expression for the atomic partition function. The operator 
P3 can be taken into account by subtracting the lowest order contributions with respect to the interaction parameter. 
As in the zero magnetic field case || one can define a Planck-Larkin partition function 

<r B (T) = (exp(-pE mQ ) - 1) + J^ 2 {exp(-0E mv ) - 1 + (3E mu ) . (35) 

Here, the factor 2 has its origin in the near-degeneracy of the hydrogen-like eigenstates. One can simplify the results 
by integrating over the pseudomomentum K 

dK z dK x exp (~^$- = (2nMkT)i^- . (36) 

Now we can rewrite Eq.(^7j). By using the eigenvalues Eo mu K z K ± ( |3l] ) and by introducing the modified fugacities 
z e j according to Eq.(Q) wc arrive at the following expression for the bound state contribution to the second virial 
coefficient 

a hound ~ ~ abound ~ ~ r. | \3 tann (gg) tanh(xj) . „-2xj\ „-2mi, M± 

m— 



M_l and cfb{T) are given by Eq.(|34|) and Eq.(|35|) ; the energy eigenvalues E mv by Eq.(|32|) and Eq.(p3), respectively. 
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V. EQUATION OF STATE 

Now we can sum up all contributions we have considered. According to Eqs. ( |l7|j20l , f23| , [37|) and expanding the ideal 
contribution in terms of the modified fugacities up to the order z 2 the pressure reads as follows 

n - k3 - - ( c ,1 7r5 tanh(a; Q ) cosh(2x a ) tt , it? , 2 

AP = 2^ Za + ^ + Z^ ZaZb I ~ SabX ab^ ^27 7 +^6 ^ KbUbfliXa) - — Kb Lb h {x a , x b ) 

a ah a 

- S ab ^ In (2) A 3 Qb e a bh(xa)) + z e ^ B*r nd ■ (38) 



(39) 



The chemical potential in Eq. (p8[) can be eliminated by using the relation 

- d((3p) 



to obtain the equation of state for a magnetized plasma. This procedure has been carried out numerically and the 
results are given in Fig. 1. 




Eqs. (^8 39) describe the ionization equilibrium of a weakly coupled hydrogen plasma in strong magnetic fields in an 
implicite form. The effect of the nonideality (i.e. of the scattering states contribution) of the plasma is to reduce the 
pressure. This contribution dominates the bound states contribution at high temperatures while at low temperatures 
the bound state term is dominant. Independent of the nonideality of the system we may also characterize the pressure 
by the magnetic field strength. For T < 6 x 10 5 K the pressure decreases with increasing magnetic field strength, while 
for T > 6 x 10 5 -ftT the pressure increases as the magnetic field increases. This can be explained by the domination 
of the lowering of the ground state energy with increasing magnetic field strength at low temperatures, while at high 
temperatures the decrease of the phase space volume dominates. 

In order to give a more explicit representation of the ionization equilibrium we will derive a generalized Saha 
equation in the next section. 



VI. SAHA EQUATION 



In previous treatments of this problem ]3|, [12[|13| ] the interaction between the charged particles has been neglected. 
But at high densities considered here interactions between the particles play an important role. Our method is based 
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on the chemical picture in which bound states are considered as composite particles, which must be treated on the 
same footing as elementary particles. By inspection of the fugacity expansion (|38| ) we reinterpret the term containing 
the partition function <r B (T) as the fugacity Zq of the neutral atoms 

2£ = ~ Zl ~z e B h ° und . (40) 

Defining the fugacities of the free composite particles in the chemical picture by z* = z e , z* = z% the pressure reads 
as follows 

f3p = z* + zt + ^-+J2 KKBlr + ~z* , (41) 

ab 

with B f a ' b ee = B s a l att + Bf b eal and al is given by 

^deai x X 3 7ri tanh(z a ) cosh(2a; a ) 

B ab -- d abA afc — -j- -. (42) 



The particle densities of the new species are given by 



*_-*9(Pp) *_~*d{[ip) *_-*»(/%>) , A ~s 

n e - z e — , n t - z t — , n -z — . [4d) 



Solving this equation by iteration we find 

\nz e =\nn e - --^—-2n e {B J ee + B J ei ), 
Inz, = Inn] - ~~ ^ + ^ ' ' 

In Zq = In rig , (44) 

where now k* 2 = (n* + n*)f3e 2 /ea = 2n*[3e 2 /ea. By inserting the fugacities according to Eq.(f44|) into Eq.(|4(]) the 
following Saha equation is obtained 

U - Ubound ( _ fc 2 ^ _ 2 * B free \ 



where B^J h ee is to be taken from Eq. ( |26| ) and (]4^) . It is useful to extend the range of validity of Eq. ( [45| ) for large £ ab 
by a kind of Pade approximation. Noting that 

2nz / B , = (1 — k a) w , (46 

ab ' 

where a may interpreted as an effective radius of the charged particles and is defined by 
Aire 2 ( i iri ir 

ab \ 

. n 7r^ tanh(x a ) cosh(2x a ) \ x/tt \— s > „ . . , „ , . . , 

+ <U X 3 ab - — A. ~ 77T E A <* b (/a^"' Xb ) + ln ( 2 )/3(^a)) , (47) 

4 £a cosn (x a J / id Qb 



we find the modified Saha equation 



The Eq. (pq) differs from the Saha equation given in || by an additional exponential factor which may be interpreted 
as the lowering of the ionization energy. In Fig. 2 the degree of ionization a = n*/n for a dense hydrogen plasma at 
various magnetic field strengths is plotted and compared with the results of the ideal Saha equation S . We find an 
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increase of the ionization degree in comparison with the ideal Saha equation [|3[ due to the nonideality effects. For 
densities of about 10 29 — 10 30 m 3 the deviation from the ideal Saha equation may be as large as 10-15% (see Fig. 2). 
At even higher densities, i.e., n 3> 10 30 to~ 3 , this result may only be used as a rough approximation. The plasma can 
no longer be regarded as a weakly coupled system rather it must be treated as a strongly coupled system. 

Additionally, we may characterize the dependence of the ionization degree on the magnetic held strength. With 
increasing magnetic field strength the ionization degree decreases at temperatures T < 6x 10 5 K, while for temperatures 
T > 6 x 10 5 K the ionization degree increases. The explanation of this effect was given in section V. 



1 




100000 1000000 
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FIG. 2. Degree of ionization at a density of p ~ 2gcm~ 3 (n = 1 X 10 30 m -3 ) for various magnetic field strength. The 
ionization fraction for k — is included (solid line). 

VII. CONCLUSION 

In this paper we constructed a theory describing a hydrogen plasma in a constant uniform magnetic field. Starting 
from a fugacity expansion we derived a general expression for the second virial coefficient as a perturbation expansion 
with respect to the interaction parameter e 2 and we explicitly calculated the lowest order contributions for the 
scattering part and considered bound state contributions at arbitrary order by using the approximate results for the 
binding energy of Lai&Salpeter M. The results were used to establish the equation of state. Finally, we have derived 
a generalized Saha equation and we have shown that at high densities and at temperatures typical for the surface of 
neutron stars nonideality effects can significantly increase the degree of ionization. 

The accuracy of the absolute values of the considered physical quantities can be improved by using more accurate 
energy eigenvalues, i.e. better fitting formulas, and by calculating even higher order contributions to the scattering 
part of the second virial coefficient. Nevertheless, the influence of the nonideality effects on the ionization equilibrium 
as shown in this paper remains approximately the same. 
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APPENDIX A: HARTREE-FOCK TERM 

By using the representation of the Green's function in terms of the spectral function (Eq.([l5])) we obtain for the 
Hartree-Fock term 
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(V )HF = \ £ Tr (CT) J^J^ £o|p ;: p2 p / /o(-0 / /oM / cffx / ^ 

x e^ lTl e^ T2 <( Pl , T x ) <(p 2j T 2 ) . (Al) 

In order to fulfill the periodicity condition of the Green's function every time variable must be extended in the 
complex time region. Therefore we associate with each time variable a small negative imaginary part t — > Hi — i8) 
and the corresponding integration may be taken in the sense of an inverse Laplace transform. Inserting Eq. ( |lq ) , the 
Hartree-Fock pressure should be written as, 

dpi f dp 2 e 2 a f , . . , f , , , . f 5+io ° ds 



2^-7 (2^)3 7 (2tt) 3 e I pi - P2 | 2 7 , ~" u ^' < / Ju ^./ a _ iDO 2™ 



PL 



^ , 5+100 2C ° Sh ^ (5l + ' S2) ) exp (-^ , , M , 1 , 

/ 5 _ ioo 2m cosh (^f Sl ) cosh (^fs 2 ) V 2 ™a / V 2m a 



This integral may be simplified in the nondegenerate case, fo(ui) — > e^e - * 3 ", where the w and s integration are 
Laplace transform and inverse, so that 



/T/\ 2 /" d Pi /" rf P2 ej cosh (cjg/3) 



2m„ 



( pIx+pIv+pIx+pIv , ,(^c Q \\ , ... 

x exp 2 a tanh -^/? . (A3) 



WaW? V 2 

The Gaussian momentum integrations are readily carried out, with the result 



nr\ ~2 \2 tanh(x a ) tt2 /" rf Pl e 2 a cosh(2x a ) / p 2 

{V)hf = > z a \ aa -3- / — — - — : - — -j- — - exp — - 

^ ia 25 j {^Y e I Pi r cosIr(a: a ) V 2 

(^2 _|_ ^2 \ 
\x V tanh ( a:a ) 1 ' ( A4 ) 

The remaining integrals with respect to P i may be evaluated exactly and the result can be expressed in terms of 
elementary functions p4[ 



(v) HF = £ - * al ^- cosh f a) -W 1 "^ . (A5) 

Finally, the charging integral may be carried out to obtain the Hartree-Fock contribution given in Eq. (|FJ 

APPENDIX B: MONTROLL-WARD TERM 

According to Eq.(|l9|), the Montroll-Ward term may be written as, 

(v)m W = i5>(^> / ^ fj^f^sf* dtv^vM 



ab 



xC(P - f ; *) G^< (p + I ; -*) (k - I ; t) Gl< (k + I ; -i) . (Bl) 

We are interested in the low density region, i.e. fo(u>) < 1. Thus we consider only contributions up to the order z 2 . 
Applying the same arguments as discussed in the previous section leads to the equation 
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(V)mw = \ E z a z b Tr M J^-j^-j^- j^A K s b (q) K*(«0 (B2) 
x A-(p - |,t) A*(p + |, -i/3 - t) Ag'(k - |, t) Ag'(k + |, - t) . 



Again, A(k) may be replaced according to Eq.([15|) and all Gaussian integrals may be evaluated with the result 

(V)mw = f E TT% ( — ) 2 A " b /' * / dq ^Tnr i exp *( J ~ f )) exp ^ - *) (v« + ' ( B3 ) 

ah ^ ' \ £ / Jo J 1 T « ^ a fc Q 

where we have defined y aj & = A^ a b& sinh(x ai f,i) sinh(x 0j b(l — t)) / (\ 2 ab t(l — 2x ai0 sinh(x a , )). Introducing spherical 
coordinates one can integrate with respect to q. The result is readily seen to be 

(V)MW = f E ^ (if) 2 / d * 1 /' ^ (^^(i - *) (»» + 1/6" ^(ife + y b - 1))) 

x (l - er/ (/€A abv /t(l - t)(y a + y& - * 2 (y a + 2/6 - 1)))) ■ (B4) 

Finally, the z- integration yields, 

A/\ P^zJk (e a e b V f 1 n 2 ( 2 exp ( K 2 Ay (1 - t) (y a + y 6 )) / r- — \ 

2 ^ ^ \ e o J Jo K V KA abx /i(l-i)((y a + ? yb-l)) V / 

2fc+l 



^E^ (2fcTT)i! — •«(?-*-5'f 1 -tT i JJ- < B5 » 

For a low density plasma we may expand this expression in powers of k\ and retain only contributions to first order. 
Using the representation of the hypergeometric function 

fl 1 3 ,\ if i arcsin(a;)\ , , 

2 F, -,^=-(yi—^ + — Al) (B6) 

the Montroll-Ward contribution to the second virial coefficient becomes 



nr\ *ST n2 ~ \ ( 1 , 4 f m /7n 7\ / _l arctanhy/1 - {y a + y b ) 

\ V )mw = kT—-^ — kTz a z b \ ab I — ~U — y dt\/t{l-t){y a + y b ) 



\/i - (y a + y b ) 



(B7) 

After performing the charging procedure one may obtain the Montroll-Ward contribution to the pressure (Eq.(EG)). 



APPENDIX C: SECOND ORDER EXCHANGE TERM 



This contribution is found to be 

<">- -5 E ^ C^I<0?I<0fJw? v ^ v ' (k) G " >(p + -2 + 

xGT (p - § + § ; -*) GTfo - § - §;*) GTfo + f - |i -*) . (CI) 

where the screened potential V s was replaced by the bare Coulomb potential V. Performing the Laplace transform 
and inverse, this equation may be rewritten as, 

= I E -2 jf ^SwrSw?Swr y(q) F(k) A * (p + f + ^ 

x<(p - f + 1; -0 - *MZ(p - f - |;*) ^(p + ~ - \\ -0 - 1) ■ (C2) 



n 



Carrying out all elementary integrals we obtain the result 

(VU =kTj2 X aa ( ) ll f 3 ( Xa ) , (C3) 



2 V47T6, 

where fs(x a ) is given by the integral representation 

1 cosh(2x a ) f f°° 1 arctanhy/v a 1 

J3{x a ) = : Zy: ttt — 7 / "* / "*i 



7rln(2) cosh 2 (x a ) Jo Jo y/h + At(l - t) Vv a hx a /[tanh(x a t) + tanh(x a (l - t))} + 1 ' 

(C4) 

with 

tx[ta,nh(x a t) + tanh(x a (l — t))]/x a + 4[tanh(a; a t) tanh(x a (l — t))]/x^ 
Va= t 1+ 4t(l-t) 

x (C5 ) 

h + [tanh (x a t) + tanh (x a (l - t))]/x a 

The charging procedure yields an additional factor 1/2 and, finally, one obtains for e 4 -exchange term the result given 
Eq.©. 
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